function [error_rate] = BER(vec1,vec2)
mismatch = find(vec1 ~= vec2);
error_rate = numel(mismatch)/numel(vec1);